Between a metal and an insulator: the critical state of the Anderson transition 
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Using a three-frequency one-dimensional kicked rotor experimentally realized with a cold atomic 
gas, we study the transport properties at the critical point of the metal-insulator Anderson transition. 
We accurately measure the time-evolution of an initially localized wavepacket and show that it 
displays at the critical point a scaling invariance characteristic of this second-order phase transition. 
The shape of the momentum distribution at the critical point is found to be in excellent agreement 
with the analytical form deduced from self-consistent theory of localization. 
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Different phase transitions observed in various fields 
of physics often share similar characteristics [l|. Of spe- 
cial interest is the behavior of the system at the critical 
point (for example scale invariance) and in its immedi- 
ate vicinity (e.g. divergence of a characteristic length 
scale). The advent of cold atom physics has offered new 
possibilities of direct experimental observation of such 
characteristics of quantum phase transitions. In this let- 
ter, we show that the Anderson metal- insulator transi- 
tion (which has only recently been observed with atomic 
matter waves [2|) obeys scale invariance at the threshold, 
defining a new state of matter between a metal and an 
insulator. 

The Anderson transition takes place in 3-dimensional 
(3D) disordered non-interacting systems in the meso- 
scopic regime (where the transport is coherent). It in- 
volves a metallic phase at low disorder associated with an 
essentially diffusive transport, and an insulating phase at 
large disorder where transport over long distance is inhib- 
ited by interference effects: this is the so-called Ander- 
son localization phenomenon [3|. The Anderson transi- 
tion is a second-order (continuous) phase transition: On 
the insulating side, the localization characteristic length 
i diverges algebraically, £ (x \K — Kc\~'^ when K, the 
control parameter, approaches the threshold K^ of the 
transition. On the metallic side, similarly, the diffusion 
constant vanishes algebraically D (x \K — Kc \ '' ■ The crit- 
ical exponents s and v are equal in 3D. and universal 
(they do not depend on the microscopic details of the 
system) [J|. Only recently have these theoretical pre- 
dictions been confirmed experimentally and the value of 
V ^ s unambiguously determined [2,l5|, l6|: z^=1.4±0.3 
is found perfectly compatible with v = 1.57 ± 0.02 ob- 
tained from numerical simulations of the 3D Anderson 
model 0,[i|. 

The state of a disordered system is, in this context, 
characterized by its transport properties. One can con- 
sider the behavior at large distances and long times of 
the (disorder) averaged intensity Green function (AIGF) 
which determines the probability P{r,r';t) for a parti- 



cle to go from r to r' in time < |i9|]. In the insulating 
phase, the AIGF is a stationary function exponentially 
localized: 



P{r,r';t) ^ exp{-\r - r'\/2e), 



[localized] (1) 



while in the metallic regime, it is a Gaussian expanding 
diffusively: 



P{r, r'; t) - cxp[-(r - r'fl2Dt]. 



[diffusive] (2) 



These two behaviors are however long time asymptotics. 
Indeed, a localized AIGF is observed only for times 
t ^ te^, where tg is the localization time (the time-scale 
associated to localization). At the transition, ti ^ i^ di- 
verges and the system becomes scale invariant. What is 
the AIGF behavior at the critical point? In the following, 
we show that it scales as: 



P{r,r';t) ^ exp 



-a\r 



_ r'|3/2/il/2 



[critical] (3) 



where a is a known (measurable) quantity. This defines 
a new state, since the shape does not change with time. 
Such a state of matter, intermediate between an insula- 
tor and a metal at all scales, has never been directly ob- 
served experimentally, although interesting results have 
recently been pub lished for ultrasound waves in the lo- 
calized regime [10| . The purpose of this letter is to report 
the first experimental characterization of such a critical 
state of the Anderson transition. 

In cold atomic gases, it is possible to prepare the sys- 
tem in a localized state and follow its evolution over 
time; this constitutes an experimental measurement of 
the (A) IGF [2, lUl, ll2| , which is impossible to achieve in 
state of art solid state physics. Observing the 3D An- 
derson transition in configuration space with cold atoms 
requires a disordered potential with a correlation length 
comparable to the de Broglie wavelength [l3| which has 
not yet been achieved. We have recently shown [2|,|5| that 
it is nevertheless possible to observe Anderson localiza- 
tion and the Anderson transition in momentum space by 
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Figure 1: First row: Measured AIGF of tiie quasiperiodic atomic kicked rotor at difTerent times t in the (left to right) localized 
{K,e) = (6,0.1), critical {K = Kc,e) = (8,0.38) and diffusive {K,e) — (11,0.8) regimes. Second row: Appropriate rescalings 
of the momentum by i° (localized), i'^" (critical) or t^" (diffusive), bring the curves at different times into coincidence (the 
vertical scales are also rescaled in order to preserve normalization). The shapes are different in the three regimes: exponential 
localization, Eq. ([l]), in the localized regime, Gaussian shape, Eq. ((2)1, in the diffusive regime, and the new "Airy shape", Eq. ((Gjl, 
at the critical point. Parameters are: S = 2.89, UJ2 = 2TTy/7, ujs — 27rVTf. Time is measured in number of kicks and momentum 
in units of 2hkL ■ 



using a different system, the atomic kicked rotor (de- 
scribed below) , where the chaotic nature of the classical 
motion replaces the disordered potential. 

Our atom-optics system (see [5| for a detailed descrip- 
tion) consists in a cloud of laser-cooled cesium atoms 
(FWHM of the momentum distribution of Shk^) inter- 
acting with a pulsed (period Ti = 27.778 ^s), far de- 
tuned (A ~ —11.3 GHz) standing wave (wavenumber 
fci = 7.4 X 10^ m~^ and one way intensity Iq ~ 150 
mW) . The amplitude of the kicks is modulated with two 
frequencies 0^2 and W3. The Hamiltonian reads: 



H 



El 
2 



JV-l 



K cos 2: [1 + £ COS (0^2^) COS (uJst)] y^ 5{t 



(4) 
where time is measured in units of Ti, space in units 
of (2fci)~^, momentum in units of 2hkL/^, with k — 
AhklTi/M = 2.89 (M is the atom mass) playing the role 
of an effective Planck constant {[x,p] = ik) and K is the 
average kick amplitude. The kicks are short enough (du- 
ration T = 0.8 /is) as compared to the atom dynamics 
so that they can be considered as Dirac delta functions. 
Decoherence processes, analyzed in detail in [5| are neg- 
ligible for the typical duration of the experiment t ~ 160 
kicks. 

If W2, W3, TT and k are incommensurate, this ID 
quasiperiodic kicked rotor has been shown to be equiva- 
lent to a 3D disordered anisotropic system [y, [ij, [l5| and 
to display an Anderson metal-insulator transition, as ev- 



idenced by the fact that it belongs to the universality 
class of the 3D Anderson model [61, |7|, i.e. has the same 
critical exponent v. Here, the localization manifests it- 
self in momentum space instead of configuration space. 
We thus expect the AIGF to take simpler forms in mo- 
mentum space, with expressions similar to Eqs. (Il])-(l3]) 
(simply replacing position r by momentum p) . In order 
to avoid confusion, we will use the notation H(p,p'; t) for 
the AIGF in momentum space. Thus, an initial momen- 
tum distribution W{p, i = 0) will on average evolve at 
time t to: 



Wip,t)^ nip,p';t)W{p',0)dp' 



(5) 



Experimentally, we are able to measure the momen- 
tum distribution at the end of a pulse sequence (up 
to 160 kicks), using Raman stimulated transitions (see 
[13, [13 for details). The initial state W{p,t = 0) is 
a thermal momentum distribution whose width is much 
smaller than the width of the final distribution, and can 
thus be approximated by a 5- function 6{p) in Eq. (0). As 
a consequence, the final momentum distribution W{p, t) 
faithfully measures the intensity Green function H(p, t) = 
H(p,0;t). 

Figure [T] shows the experimentally measured AIGF 
n(p, t) at various times in the localized, critical and dif- 
fusive regimes. The three regimes obey different scaling 
laws. In the localized regime (left column), the momen- 
tum distribution is localized - i.e. it is time-independent 
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Figure 2: Numerically simulated time-evolution of an ini- 
tially localized momentum distribution (log scale), for the 
quasiperiodic kicked rotor at the critical point of the Ander- 
son transition. The spreading follows an anomalous diffusion, 
with (p'^(f)) oc i'^", and the shape is preserved (i.e. scale in- 
variant), being neither exponentially (as it is in the localized 
regime), nor Gaussian (as it is in the diffusive regime). The 
analytic prediction, Eq. (|B|) is shown as the red thick curve 
at 1 million kicks. The agreement is excellent, without any 
adjustable parameter. Parameters are those of fig. [T] Time is 
measured in millions of kicks and momentum in units of two 
recoil momenta 2hkL- 



- at long times and thus scales as t°. In the diffusive 
regime, the average kinetic energy {p^it)) increases lin- 
early with time, so that the typical momentum scales as 
t^/^; this is manifest in the broadening of the distribu- 
tion with time seen in the right column. At the critical 
point of the Anderson transition, we observe [3, [a , as 
predicted by the one parameter scaling theory [18|, ll9| , 
an anomalous diffusion {p'^){t) ^ t^l"^ . This implies that 
the typical momentum scales as t^l^ leading to a slower 
broadening of the distribution (middle column). If the 
raw experimental data are rescaled according to these 
laws (lower row in Fig.[T|), i.e. plotted vs. pt°, -pt'^l^ and 
ptr^^^ in the localized, critical and diffusive regimes re- 
spectively, curves taken at various times coincide, which 
constitutes an experimental proof of the validity of the 
scaling laws [2J| . The shapes of the distributions are dif- 
ferent in the various regimes: exponential shape in the 
localized regime, Gaussian shape in the diffusive regime. 
The intermediate shape at the critical point is discussed 
below. 

Figure [1] is a clear manifestation of the scale invari- 
ance at the critical point. The anomalous diffusion is 
not a transient behavior and the AIGF keeps the same 
shape at the critical point. However, slightly off the crit- 
ical point, the AIGF tends gradually to either a localized 
or diffusive behavior, following the anomalous diffusion 
only for short times. To confirm this observation of scale 



Figure 3: (color online) (a) Experimental data for the rescaled 
critical AIGF (see fig. [T| averaged over time (black circles 
with error bars) and a fit given by Eq. ((6]), with p the only 
fitting parameter. The agreement is clearly excellent. The 
fitted value p is found compatible with p — r(2/3)Ac/3. The 
residual does not significantly differ from zero [panel (b)]. Fits 
by an exponentially localized (c) or a Gaussian (d) distribu- 
tion show significant deviations. 



invariance over a time scale larger than 160 kicks, we per- 
formed numerical simulations of the critical dynamics up 
to t = 10^ kicks. The result is shown in Fig. EI The ad- 
vantage of numerical simulations is that it is possible to 
explore the tails of the momentum distributions (hidden 
by noise in a real experiment). The anomalous diffusion 
- with the characteristic sub-diffusive t^^^ scaling - is 
clearly visible. Obviously, the distribution is neither ex- 
ponentially shaped (which would result in straight lines 
in the logarithmic plot), nor has a Gaussian shape (a 
parabola in the logarithmic plot). 

The form of the critical AIGF can be deduced from 
the self-consistent theory of localization [9j . This mean- 
field theory describes quantum transport in disordered 
systems at large distances and for long times |20|. It has 
been shown to be relevant for the ID periodically kicked 
rotor 



2lj and correctly predicts a metal-insulator tran- 



sition in three dimensions and the anomalous diffusion 
at the threshold: D(uA ^ (— iw)^/'^ with w the frequency 
conjugated to time [23. (the 1/3 exponent is the counter- 
part in the frequency domain of the anomalous diffusion 
(p^)(i) ^ t^^^ in the time domain). Using this critical be- 
havior, we can compute the AIGF for the quasiperiodic 
kicked rotor [23| . The details of the calculation will be 
published elsewhere; we obtain: 



n(p,i) = ^(3p3/^tr'^'Ai 



, \ -1/3 

3p'^'t] \p\ 



(6) 



where p is a parameter directly related to the criti- 
cal quantity Ac = lim (p^)/i^/^ (see P, Hi) via p = 

r(2/3)Ac/3, where F is the Gamma function and Ai(a;) is 



the Airy function. This expression is used for the plot in 
Fig. [21 The asymptotic form Eq. ([3]) comes simply from 
the limiting behavior of the Airy function for large x and 
is found perfectly intermediate between the exponential 
(localized) and the Gaussian (diffusive) shapes. 

The analytic prediction, Eq. ^, matches very well 
the shape obtained from numerical simulations of the 
quasiperiodically kicked rotor shown in Fig. [2l The only 
noticeable difference is near p = 0, where the result of the 
numerical simulation is slightly larger than the analytic 
prediction. Note that this is only observed at very long 
times, beyond 1000 kicks; this phenomenon is currently 
under study. On the time scale of the experiment (160 
kicks), this effect is invisible. 

Figure[3]shows the comparison between the experimen- 
tally measured critical AIGF and the analytic prediction, 
Eq. ^. The only fitting parameter is the global scale 
p (found in excellent agreement with p = r(2/3)Ac/3). 
Although it is visually not obvious to distinguish the ob- 
served shape from either an exponential shape or a Gaus- 
sian shape, a careful fitting procedure gives a clear cut 
result. The residual between the observed distribution 
and the analytic prediction ([6]), shown in panel (b), is 
consistently zero (within the error bars), while a fit with 
an exponential shape, panel (c), of a Gaussian shape, 
panel (d), displays significant deviations. This is fully 
confirmed by a quantitative check of the quality of the 
fit. The fit by the Airy function gives a x^ per degree 
of freedom equal to 1.09 - i.e. perfectly acceptable - 
while the exponential fit gives 4.5 per degree of freedom 
and the Gaussian fit 8.8, two unacceptably large values. 
This clearly shows that the self-consistent theory of lo- 
calization accounts for the critical AIGF and its scale 
invariance. 

In conclusion, we have studied experimentally the 
transport at the threshold of the Anderson transition. 
It obeys scale invariance, one fundamental property of 
this second-order phase transition, and this defines a new 
state, between an insulator and a metal. Its analytic 
form can be deduced from the self-consistent theory of 
localization. Work is in progress to allow experimental 
observations at longer times, which should allow us to 
characterize the small deviations observed numerically, 
whose origin could be multifractality. 
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